Cyclotron parametric resonance in microwave-irradiated two-dimensional electron 
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In recent experimental studies of the microwave-irradiated ultraclean two-dimensional electron 
systems (2DESs) in magnetic field a puzzling colossal, narrow photoresistivity spike has been dis- 
covered in the vicinity of the 2nd cyclotron resonance harmonic 2cjc. We suggest an explanation 
of this phenomenon in terms of the cyclotron parametric resonance in 2DES. We develop the hy- 
drodynamic theory of 2D electron liquid in inhomogeneous microwave electric field with frequency 
Q. We demonstrate that excitation of virtual Bernstein magnetoplasma modes at « 2uJc can 
strongly enhance an effective microwave field. That is way the fundamental mode of the parametric 
resonance is excited at very low microwave power. The resulting instability leads to the heating of 
2DES and the photoresistivity spike. 

PACS numbers; 73.43.Qt, 73.21. Fg, 73.40.-c 
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Unexpected features in the microwave (MW) response 
of the high-mobility 2DESs in the magnetic field B 
have generated a lot of interest among condensed-matter 
physics community, for a review see [l| . Transport prop- 
erties of these structures vary dramatically even at low 
MW power particular, the photoresistivity as a 

function of inverse magnetic field strongly oscillates with 
the period determined by the ratio of the MW frequency 
57 to the cyclotron frequency uJc — eB/mc where m is 
electron effective mass. With increasing MW power these 
oscillations become more pronounced with photoresistiv- 
ity at minima decreasing down to zero values 0, [H| ■ 

Recently a newphenomenon was discovered in the ul- 
traclean 2DESs [6|-|9[. Colossal narrow photoresistivity 
peak was observed at MW frequency V, close to 2uJc- 
One should also point out that all 2DESs featuring this 
peak exhibit a temperature-sensitive giant negative mag- 
netoresistance (GNMR) in the absence of the MW radi- 
ation d, H, [lol. An attempt to explain this spike was 
made in Ref. 11], but its origin still remains unclear. 

In this Letter, we propose a new approach to explain 
this phenomenon. Our approach developed within the 
hydrodynamic model of the 2DES is based on the para- 
metric resonance (PR) phenomenon in the 2DES subject 
to inhomogeneous MW pumping. The PR in the 2DES 
can be understood using analogy with a simple pendulum 
whose length varies periodically with frequency Q. The 
fundamental PR mode develops at fl equal to the dou- 
ble eigenfrequency of the pendulum. In 2DES, uic acts 
as the eigenfrequency and the fundamental PR mode is 
excited at $7 ~ 2ciJc. In the PR regime, the amplitude of 
the electron velocity increases exponentially in time and 
2DES becomes unstable. We show that this parametric 



instability takes place in the 2DES with sufficiently high 
mobility. 

The electric field acting on electrons in 2DES can 
be evaluated within the linear response theory. At fi- 
nite qRc, nonlocal effects cause an anticrossing of classi- 
cal magnetoplasmon mode and cyclotron resonance har- 
monics. Here, q is the magnetoplasmon wave vector, 
Rc ~ vf/i-^c is the electron cyclotron radius, vp is the 
Fermi velocity. This anticrossing splits the magnetoplas- 
mon spectrum into the multiple the so-called Bernstein 
modes measured in 2DESs in Refs. 12l-ll7l|. Important 



feature of these modes is opening of the narrow gaps 
(the 'Bernstein gaps') in the plasma spectrum near NlOc, 
N = 2,3, .... lifl falls in these gaps, inhomogeneous MW 
field can be strongly enhanced due to excitation of the 
virtual Bernstein modes. 

The instability leads to the heating of 2DES. The heat- 
ing destroys the GNMR, which samples exhibit in dark, 
and this results in the resistivity spike. It is worth to 
mention that PRs in 2DESs under different conditions 
were studied theoretically in Refs. pJ-21]. 

We consider microwave-irradiated 2DES positioned in 
plane z = and placed into the perpendicular magnetic 
field B = (0,0, B). We neglect the magnetic component 
of MW field and assume that the electric field is inho- 
mogeneous on the cyclotron radius scale because of the 
metal contacts to the 2DES, which significantly modify 
MW field [13]. Electric field of MW pumping is defined 
as Eo{r, t) = EQ{r) cos fit, with an amplitude En{r) de- 
pendent on the coordinates r = {x, y). For simplicity we 
consider electric field Eq^t) directed along x-axis and 
assume that it depends only on the x coordinate. 

The 2D electron liquid dynamics is described by the 
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Euler equation for the hydrodynamic velocity V = 

v{x,ty. 

dtV + - + {V,V)V = -Etot{x,t) + —[V,B], (1) 

T TO TOC 

where r is the phenomenological relaxation time, and 
an effective field Etot{x,t) is the pumping electric field 
Eo{x,t) screened (or anti-screened) by magnetoplas- 
mons. We will discuss explicit form of this field later. 
We are looking for solutions of Eq. ([ij for the high- 
mobility 2DES when fir ^ 1, lUcT ^ 1 in the electric 
field Etotix, t) = Etotix, ^) cos fit. 

The nonlinear term (V, V)V in ([T]) plays central role 
in our approach. It can be interpreted as a nonlinear, 
local and instantaneous Doppler shift of the frequency of 
excitations described by Eq. (HJ. 

In the linear approximation, solutions of Eq. ([1]) for the 
forced oscillations of velocity Vb — {Vox, Voy) can be writ- 
ten as Vox {x,t) = Vsx {x) sin fit, Voy {x,t) = Vcy (x) cos ilt, 
where 



v.Ax)^'-^^^i^,v.yix) = v.xi.f' 
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(2) 



To find a nonlinear correction dV{x,t) to the velocity 
we substitute V{x,t) = Vo{x,t) +SV{x,t) into Eq. ([T]) 
and derive the following exact system of equations for 
SV ^{SVx,SVy): 



( (dt + 1/t)SVx + V^xSVx + (Vox + SVx)SV^ 



-UJcSVy = 



'Vo.V^x 



{dt + l/T)SVy + V^ySVx + {Vox + SVx)SV^ 
+ujc5Vx = -VoxV^y 



(3) 



where prime denotes derivative with respect to x. 

The third and the forth terms in the left-hand side of 
Eqs. ([3]) as well as the right-hand side terms in these 
equations stem from the nonlinear term (V,V)V in 
Eq. (H)). The coefficients in the third and the forth terms 
in the left-hand side of Eqs. ^ periodically depend on 
time and cause the PR in the 2DES. 

The variable SVy{x,t) can be excluded from Eqs. ^ 
and the obtained equation is linearized with respect to 
SVx- We do not present an explicit form of this cumber- 
some equation. 

Following the standard procedure of the PR theory 
solution for the fundamental PR mode at fi w 2ujc can 
be written in the two-wave approximation as: 



SVx 



A{x) cos — + B[x) sm — 



(4) 



where so is the amplification coefficient, so ^ fi. The 
instability occurs at sq > 0. 

We substitute Eq. (U) into the linearized equation for 
dVx{x,t) to obtain the system of the linear equations for 



the coefficients A{x), B{x). In this derivation, we neglect 
the terms containing high order frequency harmonics as 
well as the nonlinear terms with respect to the electric 
field amplitude. 

Introducing new variables A — Vsx{x)A, B ~ Vsx{x)B 
and neglecting the difference between fi and 2ujc where 
it is possible we arrive to the following system equations: 





(5) 



where the detuning Sut = Vl/2 — and s = sq -I- 1/t. 

We exclude function B from Eqs. (O and simultane- 
ously omit the terms small by the parameter s/fi <C 1. 
For simplicity we also consider 'clean limit', assuming 
that the dissipation is sinall: s^ ^ Sijp . With these as- 
sumptions equation for A{x) takes the form: 



-di- 
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g2 y2 



A{x) = 0. 



(6) 



Eq. ([B]) is formally equivalent to the Schrodinger equa- 
tion with an effective potential energy 



U{x) ,r ^ 



y2 

* sx 



(7) 



and zero effective energy. For solutions of this equation 
to exist, the denominator in Eq. ([7]) must be sufficiently 
small while the numerator to be positive. 

Now we turn to evaluation of the screened field 
Etot{x,t) and, consequently, Vsx{x) and U{x), see 
Eqs. ([U and (O. In the linear response theory, Etot{x,t) 
is determined by the equation 
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Etot{x,t) = cos fit J 



dq ^^qx En{q) 



2ti 



e{q,ny 



(8) 



where EQ{q) is the Fourier transform of E^{x), s{q,uj) 
is the dielectric function of 2DES which includes both 
frequency and spatial dispersion. Condition e{q,uj) — 
corresponds to the excitation of magnetoplasmons with a 
dispersion law LOm,p{q)- If the function e{q,uj) is close to 
zero in sufficiently large interval of g, then virtual plas- 
mons excite and the screened electric field Etot is en- 
hanced ('anti-screening' effect). 

We use the random phase approximation to find the 
dielectric function of 2DES £{q,uj). In the collisionless 
limit, at g <C kp, fu^c ^ 27r^fcsT' ^ Ep, where hkp, Ep 
are Fermi momentum and energy, T is temperature, the 
function e{q,uj) is determined as 24-28|: 
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y^ — LtP' — iOsgnu! ' 



(9) 
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FIG. 1. The two lowest Bernstein branches Q,-{q) and ii+{q) 
of the magnetoplasmon spectrum cjmp(?). The branch fl-{q) 
has the maximum value luq at wave vector q — qo- The inset 
shows the spectrum in near the anticrossing . The width of 
the Bernstein gap is equal to 2tJc — ujo ■ 



where Vee{q) — ^ne^ / >c\q\ is the Fourier transform of the 
2D Coulomb potential, >f is a background static dielec- 
tric constant, Jn{qRc) is the nth order Bessel function of 
the first kind. Dissipation in 2DES will be included by 
introducing the phenomenological line broadening in the 
magnetoplasmon dispersion law. 

The magnetoplasmon dispersion law is determined by 
the equation e{q,io) — 0. Two lowest Bernstein branches 
of the dispersion law are shown in FiglTJ It is important 
to point out the opening of the Bernstein gap between 
Wo and 2ujc- 

We consider the situation when the frequency fl is 
positioned in this gap and the pumping electric field is 
enhanced due to the excitation of the virtual Bernstein 
magnetoplasma modes. This effect is described by the 
factor in Eq. ([8|), see also Figl2l 

We label the second and the first branches in the mag- 
netoplasmon dispersion law as ^±{q) respectively. To 
evaluate the screened electric field at the frequency 51 
close to 2lUc, we retain only the terms with n = 1,2 in 
Eq. p. 

At the next step, we estimate the wave vector qo 
corresponding to the maximum of the lower branch, 
luq = 57_((7o) and the gap Sujgap — 2a;c ~ loq. In 
the lowest order in the small parameter qRc we obtain 
qoRc = 4:.5aB / Rc, where ub = fi^x/me^ is the effective 
Bohr radius. For typical parameters of GaAs quantum 
wells = 7, TO = 0.067mo, = 3 x 10^^ cm"^) and 
LxJc/^TT — 10"'^^ Hz we obtain as ~ 5.5 nm, Rc « 0.37 
/im and q^Rc ~ 0.067 <C 1 so that the parameter q^Rc 
is really small. In the same approximation the gap is 
estimated as: 



^<^gap — 2uJc 
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Now we can estimate the screened electric field 
Etot{x^t). In this estimate we include only the resonant 
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FIG. 2. Inverse dielectric function e~^{q,Q,) describing the 
enhancement of the pumping electric field due to the excita- 
tion of the virtual Bernstein modes. The pumping frequency 
Q, lies in the Bernstein gap: 2ujc > Q. > uoq. Different lines cor- 
respond to the different values of the detuning 5uj = Q./2 — ujc- 



Suj/lOc = —1.25 X 10 (solid line), iJcj/wc 



-lO"'' (dashed 



line), 5u]/loc = —0.75 x 10 (dash-dotted line). Parameters 
of 2DES are given in the text. 



part of the screened electric field, El^^(x,t), connected 
with the maximum of e~^{q,iu) near go- 

The function e{q,uj) near qo can be represented as 
e{q, ft) = eo + {q- qo f/2M, where Eq = e{qo, fi). To es- 
timate the screened field El^^{x, t) we also use the model 
expression for the pumping field, Eq,{x) — (/So/Va;^ -I- P, 
where / is the characteristic length of the field inhomo- 
geneity. This model describes the electric field induced by 
the thin metal contact [2^. The Fourier transformation 
of the function Eq_{x) yields EQ,{q) = 21^90^^0(90' where 
Ko{ql) is the zeroth order modified Bessel function of the 
second kind. Using the asymptotics of Ko{ql) at go' *C 1 
we obtain the following expression for E'[^^{x,t): 



ipres 



{x,t) = (poy/2M/eo ln(1.12/go0 cos{qox)x 
exp{-\x\y/2Me^) cos(m). 



The value of E^^^ increases when $7 — wo and eq —>■ 0, 
see FiglSl However it saturates when jwo — il] becomes 
less than l/r*, where r* is the relaxation rate associ- 
ated with the smearing of the Bernstein gap. Therefore, 
the value of l/r* defines the maximal enhancement of 
E^St{x,t). 

An amplitude of E^^^ takes zero values periodically in 
x-space resulting in the singular behaviour of the effective 
potential energy in Eq. ([7|) at points where E^^^ {x, t) = 0. 
One can estimate the value of the pumping potential 
ipo when the effective potential becomes attractive near 
these points, i.e. the numerator of the expression in 
Eq. ([7]) is positive. For the 2DES parameters listed above 
and LUcT* ~ 10^ we calculate eo ~ 0.013, M « 0.015/ R^. 
We find very small value for the pumping potential: 
(ySo ~ 1 mV. Equation ^ has solution when the con- 
dition ipa > 1 mV is satisfied. It happens because the 
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denominator of the expression for U{x) in Eq. (O takes 
zero value while its numerator is positive. It leads to the 
asymptotic behaviour U (x) ^ — which corresponds 
the to so-called 'fall to the centre' in quantum mechanics 
[30| . In this case solutions of Eq. ^ always exist. 

One should point out that the instability develops 
at the distances from the contact of the order tt/qq ~ 
50Rc ~ 18 /im. This distance should be smaller than the 
sample's width. We suggest that in the narrow samples 
the anti-screening is not effective and the PR in such 
2DES is never realized. 

The parametric instability is developed in the samples 
with large width. However it can happen only in the 
high-mobility 2DES where the condition SujgapT* ^ 1 is 
satisfied. This condition defines the minimum frequency 
of the MW pumping. For 2DES with parameters listed 
above, ^ = er/m = 3 x 10^ cm^ / (V s) , and assuming 
T* « T, we estimate this minimum frequency as: 



2tt 



1 



2tt V 11.4 



VF 



1 



160 GHz. 



One can also estimate the minimum effective electron 
mobility needed for the instability to take place at typical 
frequency f7/27r = 100 GHz: ^;„„ = eT*/m « 13 x 10^ 
cmV(Vs). 

In summary, we have developed the hydrodynamic the- 
ory of the cyclotron PR in 2D electron liquid irradiated 
by the inhomogeneous MW field of frequency ft. We use 
linear response theory with nonlocal dielectric function 
e{q, uj) to find the screened microwave electric field act- 
ing on the 2D electrons. It is important that there is 
the Bernstein gap near 2ujc in the magnetoplasmon dis- 
persion law. If falls in this frequency gap and the 
gap width is larger than the effective relaxation rate, the 
screened electric field is enhanced due to small value of 
e(fi,q) in the gap. It is this effect that causes the para- 
metric instability to occur in sufficiently wide samples at 
the very low pumping. The PR is developed at frequency 
D./2ujc ~ 1 - b.7a%/Rl and leads to the 2DES heat- 
ing which, in turn, destroys the temperature-dependent 
GNMR observed in these samples without irradiation. 
Therefore the narrow resistivity spike arises. This mech- 
anism can explain the results of the recent experiments 
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